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ABSTRACT 

We study the propagation, reflection, and turbulent dissipation of Alfven waves in coronal holes and the solar 
wind. We start with the Heinemann-Olbert equations, which describe non-compressive magnetohydrodynamic 
fluctuations in an inhomogeneous medium with a background flow parallel to the background magnetic field. 
Following the approach of Dmitruk et al., we model the nonlinear terms in these equations using a simple 
phenomenology for the cascade and dissipation of wave energy, and assume that there is much more energy in 
waves propagating away from the Sun than waves propagating towards the Sun. We then solve the equations 
analytically for waves with periods of hours and longer to obtain expressions for the wave amplitudes and 
turbulent heating rate as a function of heliocentric distance. We also develop a second approximate model that 
includes waves with periods of roughly one minute to one hour, which undergo less reflection than the longer- 
period waves, and compare our models to observations. Our models generalize the phenomenological model 
of Dmitruk et al. by accounting for the solar wind velocity, so that the turbulent heating rate can be evaluated 
from the coronal base out past the Alfven critical point - that is, throughout the region in which most of the 
heating and acceleration occurs. The simple analytical expressions that we obtain can be used to incorporate 
Alfven-wave reflection and turbulent heating into fluid models of the solar wind. 



Subject headings: 



1. INTRODUCTION 

Beginning with the discovery of large-amplitude Alfven 
waves in the interplanetary medium four decades ago 
(Belcher, Davis, & Smith 1969; Belcher & Davis 1971), 
spacecraft measurements have demonstrated that turbulent 
fluctuations in the magnetic field and plasma velocity are 
ubiquitous in the solar wind (see, e.g., Tu & Marsch 1995; 
Goldstein, Roberts, & Matthaeus 1995; Bruno & Carbone 
2005). As in hydrodynamic turbulence, fluctuation energy 
in turbulent magnetized plasmas cascades from large scales 
to small scales. Upon reaching sufficiently small scales (typ- 
ically, of order the proton gyroradius and smaller), the en- 
ergy is dissipated, thereby heating the ambient plasma [see 
Schekochihin et al. (2009) for a recent and detailed discus- 
sion]. Turbulent heating of the solar wind is thus inevitable, 
and may play an important role in the origin and evolution of 
the solar wind (see, e.g., Coleman 1968; Barnes 1981; Holl- 
weg 1983; Tu et al. 1984; Hollweg 1986; Tu 1987, 1988; 
Hollweg & Johnson 1988; Velli et al. 1989, 1990; Grappin 
and Velli 1991; Velli 1991; Grappin et al. 1993; Marsch and 
Tu 1996; Matthaeus et al. 1999). 

The bulk of the fluctuation energy in solar-wind turbu- 
lence is at length scales greatly exceeding the proton gyrora- 
dius, is non-compressive, and to a good approximation can be 
treated within the approximation of reduced magnetohydro- 
dynamics (MHD) (Kadomtsev & Pogutse 1974; Strauss 1976; 
Zank & Matthaeus 1992; Schekochihin et al 2009). The cas- 
cade of energy to smaller scales in reduced MHD turbulence 
arises only from interactions between oppositely propagating 
Alfven waves. It has thus been known for some time that the 
launching of outward-propagating Alfven waves by the Sun 
is not sufficient to cause turbulent heating, unless some addi- 

1 Space Science Center and Department of Physics, University 
of New Hampshire, Durham, NH; benjamin.chandran@unh.edu, 
joe.hollweg@unh.edu 



tional process generates Alfven waves that propagate towards 
the Sun (in the solar-wind frame). Different mechanisms 
for producing Sunward-propagating waves have been consid- 
ered, including the parametric instability (Galeev & Oraevskii 
1963; Vinas & Goldstein 1991; Del Zanna, Velli, & Londrillo 
2001) and instabilities driven by velocity shear (Roberts et al. 
1992; Breech et al., 2008). In this paper, we focus on a third 
mechanism: the non-WKB reflection of Alfven waves. 

When the non-uniformity of the background plasma is 
taken into account, purely outward propagating waves are no 
longer a solution to the reduced MHD equations. Instead, out- 
ward and inward waves are coupled through linear terms that 
are proportional to the gradients in the Alfven speed (Velli 
1993; Hollweg & Isenberg 2007). This coupling gives rise 
to wave reflection. At sufficiently high frequencies, linear 
Alfven waves propagating in a hypothetical, steady solar wind 
are insensitive to radial variations in the Alfven speed and 
undergo little reflection. On the other hand, for wave pe- 
riods exceeding ~ 20 minutes, the transmission coefficient 
for outward-propagating Alfven waves launched from the Sun 
drops significantly below unity (Velli 1993). A number of au- 
thors (e.g., Velli et al. 1989; Matthaeus et al. 1999; Dmitruk 
et al. 2001; Dmitruk & Matthaeus 2003; Chandran et al. 
2009; Verdini et al. 2009) have discussed the propagation 
of Alfvenic fluctuations originating at the Sun, emphasizing 
the self-consistent generation of reflected waves and the re- 
sulting cascade and dissipation of wave energy. In particu- 
lar, Dmitruk et al. (2002, hereafter "D02") developed a sim- 
ple phenomenological model describing reflections and turbu- 
lence and obtained an analytical approximation for the turbu- 
lent heating rate valid close to the Sun in the limit of strong 
turbulent dissipation. Elements of their approach were incor- 
porated into detailed solar wind models by Cranmer & van 
Ballegooijen (2005, hereafter CvB05), Cranmer et al. (2007), 
and Verdini & Velli (2007, hereafter VV07). 

In view of the potential usefulness of D02's analytical ap- 
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proximation, our purpose here is to extend their approxima- 
tion out to the Alfven critical point and beyond, i.e., through- 
out the region where the principal solar wind heating and 
acceleration occur, and to obtain analytical formulas for the 
radial profiles of the turbulent heating rate and wave ampli- 
tudes. In section |2] we focus on waves with periods of hours 
and longer. In section [3] we extend our model to account for 
shorter-period waves. 



2. LOW- WAVE-FREQUENCY MODEL 

We begin with the usual equations of ideal MHD, 

dp 

dt 

3 



dt 



vV)y = -Vp b 



-V-(pv)=0, 

B VB 



4n 



- pV4>, 



and 



9B 

dt 



Vx(vxB), 
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where p, v, and B are the mass density, velocity, and magnetic 
field, 4> is the gravitational potential, p m = p +B 2 /%% is the 
total pressure, and p is the plasma pressure. We set 



and 



v = U + Sv 



B = B + 5B, 



(4) 



(5) 



neglect the Sun's rotation, and take the background flow ve- 
locity U to be aligned with the background magnetic field Bo. 
We assume that 8v and SB are perpendicular to Bo and non- 
compressive. We neglect density fluctuations 2 and take p, U, 
and Bo to be steady-state solutions of equations (fl~|i through 
(0. The vector Alfven velocity is given by 



VA 



B„ 



and the Elsasser variables are defined as 



8v =f §b, 



(6) 



(7) 



where 8b = 8B/ ^/4np. Re-writing equations (O and (f3]) in 
terms of z±, we obtain (Velli 1993; VV07) 



8z± 
~df 



(U±v a )-Vz ± +z : ' : -V(Utva) 



V-v A T2 V 'U 



•Vz d 



(8) 

The Vp m term cancels the compressible part of the remaining 
terms in equation (0 to maintain the incompressibility condi- 
tion V ■ z* = 0. 

To proceed further, we assume that the background mag- 
netic field Bo possesses a field line that is purely radial. This 
field line can be at any heliographic latitude. We then con- 
sider a cylindrical coordinate system (/f , <j>, z) whose z axis is 
aligned with this radial field line, and we restrict our analy- 
sis to the region close to the z axis, in which we assume that 

2 To make this assumption self-consistent, we must also neglect spatial 
variations in the background density in the directions perpendicular to Bo- 



= I/a = and dU = dv^/d^ — 0. Near the z axis, the 
components of U and va can be expanded in powers of R, 



U R =ai(z)R- 
VA,R = a2(z)R- 



(9) 
(10) 



Keeping just the first terms in these Taylor series, we find that 
dU R /dR = U R /R, dv A . R /dR = v A . R /R, and 



z t.V(Utv a ) 



R 



We define 



5- = (Bo/Bo) • V, 
as 



(11) 
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so that s is distance along a background magnetic field line. 
Since VR = R, 



UrTvam _ (Utva) aR{s) 



R 



R 



(13) 



where R(s) is to be interpreted as the distance from the z axis 
to a given background magnetic field line. 

We define vector versions of the Heinemann & Olbert 
(1980) variables f and g as follows: 



z 1 = 



z = 



1/4 



l+ri'/z 

fr, 1 / 4 
l-r,V2' 



where 



11 = p/p«, 



(14) 
(15) 

(16) 



and p„ is the value of p at the Alfven critical point (on the z 
axis). Mass conservation requires that d(pU /Bo) /ds — 0, and 
thus r) 1 / 2 = va/U. With the use of equations ( ITTb and ( [13] ), 
equation ^ can be re-written as (VV07) 3 



— + (U + v A ) + [U + v A )f 



Tl 



1/4 



z ■ Vz -1 



Vpto 



(17) 



df , df , , dL 

S + (^-VA) S + (tA-VA)g 5 



V /2 -i 



11 
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z+ • Vz" 



(18) 



where L — ln[R(s)r\ l / 4 / R le {], and R le f is an arbitrary constant. 
Close to the z axis, d(B Q R 2 )/ds = 0, and R lef can be chosen 
so that 



L= --In 
2 



V'A 
VAa 



(19) 



where v Aa is the Alfven speed at the Alfven critical point on 
the z axis. 

Dmitruk et al. (2002) numerically solved equation ([8]l in 
the limit U — * 0. They found that the numerical results could 

3 Scalar versions of the linear terms in equations {T7J and U8t (for just the 
(|> components of g and f) were obtained by Heinemann & Olbert (1980), who 
took the fluctuations to be axisymmetric with z (|>. 
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to a large extent be understood in terms of a phenomenologi- 
cal model in which z^-Vz* is replaced by |z =F |z ± /2A,(s). We 
employ the same phenomenological model [see Chandran et 
al (2009) for a more detailed discussion], replacing the non- 
linear terms z^ • Vz* + p~ l Vp tot in equations (fTTI i and (H~8b 
with |z ± /2A,(j'). Equations (T% and ([T8T l then become 



| + (^ + VA)| + (f/ + V A )f| = . 



1 



1/4 



2A, 



/ 



1 T|!/2 



g (20) 



and 



3f . , 9f 



Tl 1/4 |^|f 



(21) 



Observations (e.g. Belcher and Davis 1971) show that most 
of the Alfvenic power in the solar wind is at long periods, of 
the order of hours, as seen in the spacecraft frame. In this 
section, we thus restrict our attention to waves with periods of 
hours and longer in the frame of the Sun. We further assume 
that the amplitude of the outward propagating waves is suffi- 
ciently large, and A, is sufficiently small, that the time scale for 
the cascade and dissipation of z~ fluctuations, A/z + , is shorter 
than the wave period. In this case, we can neglect the first and 
second terms on the left-hand side of equation (|2T1 i and solve 
equation (|2~TT > to obtain 



f 



2X(i\- 



Tl 



1/4 



\)U dL 
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(22) 



where g — g/g. We assume that the rapid cascading of z 
energy leads to the inequality 



z <z^ 



(23) 



so that we can neglect the third term on the left-hand side of 
equation (1201 . Taking the dot product of equation d20l with g, 
and then taking the time average of the resulting equation, we 
find that 



-2(g 2 



dL 

a7 



(24) 



where (. . .) denotes a time average. Since we restrict our anal- 
ysis to the vicinity of our z axis, s ~ r, where r is the radial 
coordinate in spherical coordinates centered on the Sun. We 
take L to be a decreasing function of r near the Sun, to have 
a minimum at r = r m (corresponding to the maximum of va), 
and to increase with r at r > r m . We take r m to be less than 
the radius r u of the Alfven critical point. Integrating equa- 
tion d24"i >. we find that 



,2x1/2 : 



V'A 
VAa 



1/2 



for r>r m , where g a is the rms value of g at r = r a . For r < r„ 
we find that 

^\l/2 =p .f V Am ) (26) 



{g 2 ) 1/2 =gJ^-) , 

\VAaVAj 



where VA m is the Alfven speed at r = r m . 

The energy density of the Alfvenic fluctuations is p [(z + ) 2 + 
(z~) 2 ]/4. We assume that the energy that is drained from z + 
and z~ fluctuations by the nonlinear "damping" terms in our 
model is converted to thermal energy. The turbulent heating 
rate, for z~ <C z + , is thus 



Q = 



p\z-\{z 



+ \2 



41 



(27) 



With the use of equations 
average heating rate is 



, and ( 1261 1. we find that the 



Q 



pn 1/2 g 2 a u 

4(l+Tll/2) VAa 



dVA 



dr 



for r > r m , and 



m l/2 glu 

4(l+T| 1 /2) 



2 



dVA 

dr 



(28) 



(29) 



for r < r m . 

To check the validity of some of our approximations and 
compare to observations, we consider a model solar wind in 
which the proton number density is given by equation (4) of 
Feldman et al. (1997), which describes coronal holes out to 
several solar radii, plus an additional r~ 2 component: 



n{r) 



3.23 x 10 8 2.51 x 10 6 1.85 x 10 s 



J5.6 



,.3.76 



(30) 

where x = r/R Q and R Q = 1 solar radius. This leads to n = 
4 cm~ 3 at 1 AU. We take the magnetic field strength to be 
(Hollweg & Isenberg 2002) 



Bq = 



1.5(/max-l) 



1.5 



Gauss, 



(31) 



where / max is the usual super-radial expansion factor, which 
we set equal to 5. The solar wind speed U is given by flux 
conservation: 



U = 9.25 x 10 lz -cms 
n 



(32) 



where B is Bq in Gauss and h is n in units of cm~ 3 . Equa- 
tion fl32"l ) implies a flow speed of 750 km s _1 and a proton 
flux of 3 x 10 8 cm~ 2 s~' at 1 AU. The Alfven critical point 
in this model is at r a = 11.1/?©, and the maximum of va is 
at r,„ = 1.6O2?0. The solar wind speed and Alfven speed are 
plotted in Figure Q] 
We set 



X = Xo 



B(R e ) 



1/2 



(33) 



where B(R & ) — 7.5 G from equation d3Tl ), and A-o is a constant. 
Previous studies suggest a range of possible values for Xo- 
One choice is to set A,o equal to the mean spacing between 
photospheric flux tubes. Following Spruit (1981), we take this 
mean spacing to be R p ^/%B p /B aV g, where R p = 200 km and 
(25) Bp — 1 500 G are the typical radius and magnetic field strength 
of a photospheric flux tube, and B avg is the average magnetic 
field at the coronal base above the flux-tube merging height. 
If we set Bavg equal to be the value of Bo in equation OTT ) as 
x — * 1 (which is 7.5 G), then A-o = 5000 km. On the other hand, 
D02 used Xo = 3 x 10 4 km (comparable to the average size of 
supergranules), CvB05 considered values of A-o in the range 
1 — 30 x 10 3 km, and the quantity corresponding to our Xq 
in VV07 was 1.7 x 10 4 km. We will consider several values 
of Xq in the range spanned by the above values. 

To check whether we are justified in neglecting the terms 
that we have dropped in equations ( f20b and (l2Tb . it would 
be desirable to compare our approximate solutions to exact, 
time-dependent solutions of equations (f20b and (|2TT >. How- 
ever, such time-dependent solutions are beyond the scope 
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FIG. 1. — The sola r wi nd veloci ty ( solid line) and Alfven speed (dashed 
line) from Equations {30} through j32) with / milx = 5. 



of this study. Instead, we numerically solve the time- 
independent versions of equations (f20b and (fJTJ, taking / = 
at the Alfven critical point, setting A-o = 5000 km and g a = 
8 x 10 7 cm/s, and assuming that the directions of f and g are 
the same and do not change with s. This steady-state numeri- 
cal solution is still a useful reference point since we are con- 
sidering only long-period oscillations in this section. In Fig- 
ure |2] we plot the numerical solutions for |/| and g, as well 
as our approximate analytical solutions from equations (l22l i. 
d25l l, and ( l26l >. For r > 2J? Q , our approximate solutions are 
close to the numerical solutions. Closer to the Sun, our an- 
alytical expressions over-estimate the values of |/| and g. In 
particular, our approximate value of g is 23% larger than the 
numerical solution at r — R@. If A-o is increased to 10 4 km 
(1 .5 x 10 4 km), the approximate solution is 43% (57%) larger 
than the numerical solution at r = R & , although the approx- 
imate solutions remain reasonably accurate at r > 2Rq. In- 
creasing Xq reduces the nonlinear terms in equations (l20l and 
(fJTJ and degrades the accuracy of our approximations. 

The loss of accuracy in our analytical approximations as 
Xq increases is further illustrated in Figure [3] which is analo- 
gous to Figure 5 of D02. The dashed line in this figure shows 
Q/p from equations d28l and i29\ . The three solid lines plot 
the total heating rate Q total = p[(z+) 2 |z~| + (z~) 2 |z + |]/4A, in 
the numerical solutions of the time-independent versions of 
equations (l20b and (fJTJ for g a = 8.0x 10 7 cm/s and for three 
different values of Xq. Going from the top of the figure to the 
bottom, these values are 2500 km, 5000 km, and 10 4 km. At 
r > 2Rq, equation (|28l overestimates the numerical heating 
rate by < 50% for the chosen values of A-o- At r < 1.6Rq, 
equation (|29| i overestimates the numerical heating rate by al- 
most an order of magnitude in the worst case (Xq = 10 4 km). 

In FigureHJ we use our time-independent solutions of equa- 
tions @ and (ED with g a = 8.0 x 10 7 cm/s and A, = 5000km 
to plot each term in equation (fJTJ, divided by (1 — t] 1 ^ 2 ) to 
make the plot easier to read. Although we dropped df/dt 
to obtain these time-independent solutions, we estimate the 



FIG. 2. — A comparison between our approximate solutions for / | (do tted 
line ) and g (short-dashed line) and numerical solutions to equations )20t and 
)20 obtained in the limit of zero frequency. The solid line is the numerical 
solution for g and the long-dashed line is the numerical solution for |/|. For 
these curves, we have taken g d = 8.0 x 10 7 cm/s and ^-o = 5000 km. 
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FIG. 3. — The dashed line is the analytical approximation of the heating rate 
per unit mass Q/p from equations 128) and {29}, with g a = 8.0 x 10 7 cm/s. 
The solid lines are the total heating rate Qtotal — P[(z + ) 2 |z~l + (z~) 2 |z + |]/4A, 
divided by p in numerical solutions of the time-independent versions of equa- 
tions {20} and )2U . with g a = 8.0 x 10 7 cm/s and with three different values 
of Xo: 2500 km, 5000 km, and 10 4 km (going from the top curve to the bottom 
curve). As ^0 increases, our approximation becomes increasingly inaccurate. 



approximate magnitude of the df/dt term by multiplying our 
numerical solution for / by 2k/ P, where P is the wave period, 
which we take to be three hours. Our estimate of the df/dt 
term is smaller than the terms that we have kept inside the 
Alfven critical point, but df/dt increases relative to the dom- 
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FIG. 4. — The different terms in equation 12 IK divided by (1 — T] 1 ' 2 ) ( see 
text). The solid line is the term proportional to dL/ds, the l ong-dashed line is 
the nonlinear term on the right-hand side of equation 12 it . the dotted line is 
the d//df term, and the short-dashed line is the 3/ /ds term. As in figure[2] 
we have taken g a = 8.0 x 10 7 cm/s and Xq = 5000 km. 



FIG. 5. — The thin solid line is the value of Q in equations (28) and i29i 
for g- d = 8.0 x 10 7 cm/s . The dashed line is the total (electron plus proton) 
heating rate in model SW2 of Allen et al. (1998). The thick solid line is the 
proton heating rate from Eq. (14) of Cranmer et al. (2009), multiplied by 4/3 
to convert it to an approximate total heating rate (the proton heating rate is 
typically 60-90% of the total heating rate in their model). 



inant terms as r increases towards r a . We have compared the 
different terms in equation (l2lT > at larger r, and find that our 
neglect of the df/dt term becomes unjustified for r > 2QRq, 
because the nonlinear time scale becomes larger than the wave 
period (measured in the Sun's frame) due to the increase in A, 
and the decrease in z + as r increases. Figure |4] shows that 
the df /ds term, which we have neglected in equation (ETt . 
is negligible except in a small region near r = r m , where our 
solution for / passes through zero. Although it is not plotted, 
we have also carried out the same comparison for the different 
terms in equation d20t . using our time-independent numerical 
solution. We find that the dL/ds term that we neglect in equa- 
tion d20t is indeed smaller than the other terms out to well 
beyond 20 R & , because / <C g. 

In Figure \5\ we plot the turbulent heating rate in equa- 
tions d28l and ( |29l (thin solid line). Although the approx- 
imations leading to equations d28l and ( |29l break down at 
r > 20R & , we plot equation (|28T i out to beyond lOORp to illus- 
trate how the analytical formulas extrapolate to larger radii. 
The heating rate in our model depends on a single free pa- 
rameter, g a , which we set equal to 8.0 x 10 7 cm/s in order to 
roughly match two previous empirically constrained models 
of plasma heating in the fast solar wind (Allen et al. 1998; 
Cranmer et al. 2009). We note that for zT <C z + , this value 
of g a corresponds to a fluctuating velocity of 200 km/s at the 
Alfven critical point. The dashed line in Figure [5J is the to- 
tal (electron plus proton) heating rate £>A98 in model SW2 of 
Allen et al. (1998). The thick solid line is the proton heat- 
ing rate in the model of Cranmer et al. (2009) (their Eq. 14, 
which applies at r > 0.3 AU), multiplied by 4/3 to convert the 
proton heating rate to an approximate total heating rate, de- 
noted 2co9- [The proton heating rate is typically 60-90% of 
the total heating rate in Cranmer et al's (2009) model]. 

Model SW2 of Allen et al. (1998) does not provide a tight 
constraint on the heating rate, since widely different heat- 



ing functions and momentum-deposition profiles (their mod- 
els SW2, SW3, and SW4) were shown to be consistent with 
in situ observations of the fast wind. We include 2a98 m Fig- 
ure [5J to show that our results for the heating rate at 2Rq < 
r < 7Rq are reasonably close to at least one empirically con- 
strained model for a plausible value of g a . On the other hand, 
our heating rate greatly exceeds £?A98 at r < 1 .5Rq, This dis- 
agreement may be due in part to the effects described in sec- 
tion |4] and also to the point illustrated in figure [3] that we 
substantially overestimate the heating rate close to the Sun 
when A-o > 5000 km. We emphasize that although Figure [5] 
represents a rough first attempt to compare our model to ob- 
servational constraints on the heating rate, this comparison 
must be viewed with some caution, because the density pro- 
files in the three models are different, and Q depends upon p 
in equations (f28b and i29i . 

In Figure [6] we plot the rms amplitude of the fluctuating 
velocity in our low-wave-frequency model as a function of r 
(dashed line), using equations d25l l and d26l i and the relation 
8v 1Tns ~ ((z + ) 2 ) 1/2 /2, which is valid when z~ <§; z + - We also 
plot observed values and upper limits from remote UVCS 
observations of the corona (Esser et al. 1999; CvB05) and 
in situ measurements from Helios and Ulysses (Bavassano 
et al. 2000). (The solid line in this figure is from the ex- 
tended model described in section[3]) For g a = 8 x 10 7 cm/s, 
our low-wave-frequency model over-predicts the rms ampli- 
tude of velocity fluctuations at l.5R@ < r < 2R & . As noted 
above, the approximations leading to equations (l25l l and d26i > 
break down at r > 2QR & . Nevertheless, we plot equation d25T > 
out to larger r to illustrate how the analytical formulas ex- 
trapolate to larger radii. This extrapolation shows that for 
g a = 8 x 10 7 cm/s, equations (1231 1 and (|26| | under-predict 8v'nns 
atr >70R e . 
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FIG. 6. — The dashed line is the rms value of the wave velocity in our 
low-wave-frequency model, which is approximately z + /2 when r « z + , 
obtained from equations \25\ and (26) with g„ = 8 x 10 7 cm/s. The solid 
line is the rms value of the wave velocity in the extended model described in 
sectionfJJ obtained from equations (57) and (38) with g a = 7.2 x 10 7 cm/s and 
X = 0.65. The open box is adapted from Figure 9 of CvB05, and represents 
an upper limit on the non-thermal line widths obtained by Esser et al. (1999) 
using off-limb UVCS measurements (Esser et al. 1999). The open diamonds 
are from in situ Helios and Ulysses measurements (Bavassano et al. 2000). 



3. EXTENDED MODEL ACCOUNTING FOR SHORTER-PERIOD 
WAVES 

In section[2] we took the waves launched by the Sun to have 
periods of hours and longer. On the other hand, CvB05 sug- 
gested that a significant fraction of the outward wave power 
near the Sun is in waves with periods of 1 minute to 1 hour. 4 
In this section, we extend our model to account for such 
shorter-period waves in an approximate manner. The approx- 
imations we use are very simplistic, but are not limited to 
r<20R e . 

In the spirit of our simple, phenomenological modeling, we 
begin by replacing each g and f in equations d20l > and OTT l 
with the (time-independent) rms values (g 2 ) 1 / 2 and (/ 2 )'/ 2 , 
respectively. As shown by Velli (1993) and others, wave re- 
flection becomes less efficient as the wave period decreases 
below ~ 1 hour, since higher-frequency waves with shorter 
parallel wavelengths are less affected by the radial gradient in 
the Alfven speed. To account for this effect, we set 



dL 



dL 



ds ^ ds 



(34) 



in equations d20i > and (ETT i. where % is a dimensionless con- 
stant. In the low-wave-frequency limit of section [2J % — > 1. 

4 We note, however, that VV07 considered two different frequency spectra 
for the z + waves at the coronal base: a flat spectrum in which the energy 
is mostly in waves with periods P of ~ 10 minutes (the small-P cutoff in 
their spectrum), and a steep spectrum in which most of the energy is in the 
longest-period waves in their model, with P ~ 10 3 hours. With the steep 
spectrum, their model agreed with in situ measurements of the z + and z~ 
amplitudes, but with the flat spectrum their model over(under)-estimated the 
z + (z~) energy relative to the in situ measurements. The VV07 results thus 
favor a picture in which the Sun launches primarily long-period waves. 



To model the mix of wave frequencies launched from the Sun, 
we leave % as a free parameter between and 1. We take 
(g 2 ) 1 / 2 and (/ 2 ) 1,/2 to have scale lengths in the direction of 
Bo of order r. We continue to assume that / <C g, and set 
(r| - l)dL/ds — > |(r] — l)dL/ds\ in equation (l2lT) to maintain 
the positivity of (f 2 ), so that equation d2TT > becomes 



2\l/2 



2%X\r\-l\U 
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1/4 



ds~\ 



(35) 



Dropping the third term on the left-hand side of equation d20i > 
and using equation (|35t . we obtain 



which we integrate to obtain 



dL 
ds~ 



(S 



2x1/2 _ 



for r > r,„ and 



Or 



2\l/2_ 



VA 
VAc, 



"Am 



JC/2 



X/2 



for r < r m . From equation (1271 . we then find that 



for r > r,„, and 



Q 



4va(1 



Xm l/2 8 2 a U 





dv A 


\VAa J 


dr 



4v A (l +r| 1 /2) \v Aa V K 



Am 



(36) 



(37) 



(38) 



(39) 



(40) 



dvA 
dr 



for r < r m . The results of our low-frequency model in sec- 
tion[2]are recovered by setting % = 1 in equations ( f35l > through 

When <C z + , the rms amplitude of the fluctuating ve- 
locity Svrms is approximately ((z + ) 2 ) 1//2 /2, and there are two 
free parameters that determine Sv^: g a and %. We set 
g a = 7.2 x 10 7 cm/s and % = 0.65 in order to match the obser- 
vational constraints on the the fluctuating velocity, shown in 
Figure [6] The solid line in Figure [6] is the resulting value of 
Svjjns obtained from equations d37| i and (1381 , In Figure Q we 
plot the heating rate from equations ( |39i > and (l40l for these 
values of g a and X, along with the empirically constrained 
heating rates Qa98 and Qco9 from the models of Allen et 
al. (1998) and Cranmer et al. (2009) described in section [2] 
We plot the rms amplitudes of z + and z~ resulting from equa- 
tions ((35), d37l i, and d38l in Figure [8j as well as in situ mea- 
surements from Helios and Ulysses (Bavassano et al. 2000). 
The rms amplitude of z~ that results from equation (l33T > de- 
pends upon both % and Xq. To match the in situ data, we set 
A.o = 1.6 x 10 4 km in this figure, keeping % = 0.65 as in Fig- 
ures|6]and|7] As Figure [8] shows, the approximation z~ <C z + 
fails at r > 1 AU. 

Overall, the model matches the observations and previous 
models fairly well. However, as in our low-wave-frequency 
model, the heating rate in our extended model is much larger 
than (2a98 at r < 1 .5Rq, This discrepancy may arise in part for 
the reasons discussed in section [4] In addition, our extended 
model is likely similar to our low-frequency model in signifi- 
cantly overestimating Q close to the Sun when A-o > 5000 km 
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FIG. 7. — The solid line is the value of Q in equations (39) and fl40t for 
g a = 7.2 x 10 7 cm/s and % = 0.65. The dashed line is the total (electron plus 
proton) heating rate in model SW2 of Allen et al. (1998). The dotted line 
is the proton heating rate from Eq. (14) of Cranmer et al. (2009), which is 
typically 60-90% of the total heating rate in their model. 



in the three models are different, and Q depends upon p in 
equations ( [39b and d40b . 

Compared to our model for low-frequency waves in sec- 
tion 12 our extended model accounting for shorter-period 
waves does a notably better job of reproducing the obser- 
vational constraints on the velocity fluctuations in the solar 
wind. This is not surprising, given that there is an additional 
free parameter in the extended model and the approxima- 
tions underlying our low-wave-frequency model break down 
at r > 2Qr?0. There is, however, an additional reason for the 
closer agreement between the extended model and the data. 
As can be seen in Figure|6l the value of 8v rms in our low- wave- 
frequency model drops off more rapidly with radius than does 
the data. When the efficiency of reflection is reduced, the rate 
at which the outward wave energy cascades and dissipates is 
reduced, and so \(d /dr)8vrms\ decreases. 

4. THE ENERGY CASCADE TIME FOR Z+ WAVES 

If the Sun launches an outward-propagating Alfven 
wave packet with a perpendicular scale length A-o at the coro- 
nal base, with A-o greatly exceeding the dissipation scale, 
then this wave packet must propagate some distance into 
the corona before its energy cascades to small perpendicu- 
lar scales and dissipates. We estimate this distance using ap- 
proximations and results from our analytical models. As in 
equation d20b . the rate at which z~ waves shear z + waves is 
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FIG. 8. — The solid line is therms amplitude of z + for g„ = 7.2 x 10 cm/s 
and x = 0.65 from equations {57} and (38}. The dashed line is the rms am- 
plitude of z~ from equation {22} for Xq = 1.6 x 10 4 km and the same value 
of X- The diamonds and crosses are the rms values of z + and z~ , respectively, 
from Bavassano et al. (2000). 



(Figure [3]). A second difference between the heating rates 
plotted in Figure [7] is that Q becomes smaller than Qcm at 
r > 1 AU. This discrepancy may arise because of the increas- 
ing importance at large r of sources of inward waves (such as 
velocity shear) that are not included in our model. We note, 
however, that some caution is warranted when comparing our 
heating rate to Qa9H and 2co9, because the density profiles 



The quantity 



\dn 



2X(U + v A )' 



(41) 



(42) 



is the time integral of Y nl [r(f),f] evaluated at a position r(f) 
that moves outward with the z + wave packet at speed U + va. 
starting from the time the wave packet leaves the coronal base 
at r = rj, ~ IRq until the time the wave packet reaches ra- 
dius r. Roughly speaking, the z + wave packet must propagate 
out to a distance at which T > 1 before its energy cascades and 
dissipates. We evaluate the right-hand side of equation d42l 
in our low-frequency model using equation d22b . To evaluate 
equation d42l using our extended-model results, we replace z~ 
with ((z - ) 2 ) 1 / 2 and use equation (f35T >. The integral in equa- 
tion d42l can then be evaluated analytically to yield 

(x/2)ln(v A /v A fo) \ir <r m 



(x/2)ln[v 2 ffl /(y Afo y A )] ifr>. 



(43) 



where the result for our low-frequency model is obtained by 
setting % = 1 • The value of T is plotted in Figure [9] assum- 
ing X = 1 and using the profiles for n, Bq, and U given in 
equations (f30b through d32l . The condition T > 1 only holds 
for r > IRq, suggesting that the turbulent heating rate may 
be smaller than in our model within a few R & of the solar 
surface. We note, however, that the reduction in Q due to the 
condition T < 1 may itself be mitigated by a further complica- 
tion, that turbulence within the chromosphere [where z + ~ z~ 
(CVB05)] may lead to the launching of Alfven waves into the 
corona with a broad range of perpendicular length scales, ex- 
tending to very small values (Chandran 2008). Further work 
is needed to explore how these effects modify the radial pro- 
file of the turbulent heating rate. 
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FIG. 9. — The time integral of the shearing rate experienced by a z + wave 
packet between the time it leaves the coro nal base at r = ~ 1 A' . and the 
time it reaches radius r, from equation (43) with % = 1 . 

In one sense, Figure [9] provides a rough consistency check 
on models in which the solar wind is heated and accelerated 
by low-frequency Alfven-wave turbulence. If F were <C 1 
at r = r a , then only a tiny fraction of the Alfven-wave en- 
ergy launched from the base of the Sun would dissipate in 
the region where most of the heating and acceleration of the 
solar wind takes place. A small fraction of the Alfven-wave 
energy would be unable to generate the solar wind, because 
the Alfven wave energy flux measured in situ at r > 0.3 AU 
is much smaller than the total solar-wind energy flux, even 
allowing for loss of Alfven-wave energy due to work done 
on the flow by the ponderomotive force. On the other hand, 
if r rose from at r = ry to a value 3> 1 at say r = 7R Q , 
then most of the wave energy would dissipate very close to 
the Sun (within the radius at which F reached a value of a 
few). If the heating rate were to fall too rapidly with r, then 
the heating profile would become inconsistent with UVCS ob- 
servations that show ion temperatures staying flat or increas- 
ing with radius out to at least 3>.5Rq, despite adiabatic cooling 
(Kohl et al 1998; Antonucci, Dodero, & Giordano 2000). 

5. CONCLUSION 

In this paper, we consider Alfven-wave reflection and tur- 
bulent heating in the solar wind using the Heinemann-Olbert 
equations, which describe non-compressive MHD fluctua- 
tions in an inhomogeneous medium with a background flow 
parallel to the background magnetic field. We approximate 
the nonlinear terms in these equations using a simple phe- 
nomenological model from D02. Wave reflection plays an 
essential role in our calculation of the turbulent heating rate, 
because only outward-propagating waves are generated by 
the Sun, and turbulent interactions arise only from interac- 
tions between oppositely propagating waves (for the low- 
frequency, non-compressive Alfven waves that we consider). 

In section|2] we restrict our attention to waves with periods 
of hours and longer. Our focus on low-frequency waves in 
this section is motivated by in situ observations showing that 
most of the wave power is at periods of hours in the space- 



craft frame. We assume that the energy in inward-propagating 
waves is much less than the energy in outward-propagating 
waves, i.e., z~ <C z + , an approximation that is appropriate for 
coronal holes and the fast solar wind at r < 1 AU. We also 
assume that the wave period is longer than the energy cas- 
cade time for Sunward-propagating waves, ~ X/z + . With the 
use of these approximations, we solve our model equations 
analytically to obtain expressions for the rms amplitudes of 
outward and inward waves as a function of r. We also ob- 
tain an analytical expression for the radial profile of the tur- 
bulent heating rate. There is one free parameter in the heating 
rate in our low-wave-frequency model: the amplitude of the 
outward-propagating waves at some fixed reference point. We 
take this reference point to be the Alfven critical point r = ta, 
which is at \\.\Rq in our model solar wind. Interestingly, as 
in the model of D02, the heating rate does not depend upon 
the choice of the dominant perpendicular length scale of the 
turbulence, X. Ordinarily, the heating rate is inversely pro- 
portional to X. However, in our model, the amplitude of the 
inward-propagating waves is proportional to X, because small 
X leads to rapid turbulent dissipation of the inward waves. 
Asa result, (2~p((z + ) 2 )((z~) 2 ) 1//2 /4A, becomes independent 
of X, as in equation (5) of D02. Our model can be thought of 
as generalizing the phenomenological model of D02 by ac- 
counting for the solar wind velocity, so that the model can be 
applied all the way from the coronal base out past the Alfven 
critical point. 

As discussed in section|2] one of the approximations in our 
low-wave-frequency model — that the energy cascade time of 
Sunward-propagating waves is shorter than the wave period 
measured in the Sun's frame — breaks down for r > 20R & 
for a wave period of ~ 3 hours. In addition, our low-wave- 
frequency model does not account for waves with periods 
much shorter than one hour, which may make a significant 
contribution to the outward wave flux from the Sun (CVB05; 
but see VV07). To overcome these limitations, we develop an 
"extended" model in section [3] that approximately accounts 
for waves with shorter periods and does not require the non- 
linear time scale to be shorter than the wave period in the 
reference frame of the Sun. The main additional ingredient 
in this extended model is a new free parameter, that mod- 
els the decrease in the efficiency of wave reflection at shorter 
wave periods. For appropriate choices of the two free param- 
eters in the model (% = 0.65 and g a = 7.2 x 10 7 cm/s), our 
extended model is successful at matching observations of ve- 
locity fluctuations and the Elsasser fields. 

Both of the models we have developed provide analytical 
expressions for the turbulent heating rate that can be used to 
incorporate Alfven-wave reflection and turbulent heating into 
fluid models of the solar wind, both inside and outside the 
Alfven critical point. We conclude by summarizing the rela- 
tive advantages and disadvantages of the two models: 

1. The low-wave-frequency model described in section [2] 
involves a simplistic phenomenological approximation 
of only the nonlinear terms in the basic equations. The 
heating rate in the model involves only a single free pa- 
rameter, g a , which determines the amplitude of the out- 
ward waves at the Alfven critical point. On the other 
hand, one of the assumptions of the model breaks down 
at r > 2QR & , where the nonlinear time scale becomes 
longer than the wave period measured in the frame of 
the Sun (for wave periods of ~ 3 hours). In addition, the 
model does not account for waves with periods much 
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shorter than 1 hour, which may contribute significantly 
to the outward wave flux (CvB05). 

2. The extended model described in section[3]does a better 
job of matching observational constraints on the wave 
amplitudes, accounts (in a very approximate way) for 
both long-period waves and waves with periods much 
less than one hour, and can be applied all the way out 
to 1 AU. On the other hand, this model involves sim- 
ple phenomenological modeling of both the nonlinear 
terms and the linear reflection term. The heating rate in 



the model also involves two free parameters instead of 
one. 



We thank the anonymous referee for a very helpful re- 
port. This work was supported in part by NSF Grant 
ATM-0851005, by DOE Grant DE-FG02-07-ER46372, by 
NSF-DOE Grant AST-06 13622, and by NASA Grants 
NNX07AP65G and NNX08AH52G. 



REFERENCES 



Allen, L. A., Habbal, S. R„ & Hu, Y. Q. 1998, J. Geophys. Res., 103, 6551 
Antonucci, E., Dodero, M. A., & Giordano, S. 2000, Solar Phys., 197, 1 15 
Barnes, A. (1981), Turbulence and Dissipation in the Solar Wind, in Solar 
Wind Four, edited by H. Rosenbauer, p. 326, Max Planck Institut fr 
Aeronomie Report No. MPAE-W-100-81-31, Katlenburg-Lindau. 
Bavassano, B., Pietropaolo, E.. & Bruno, R. 2000, J. Geophys. Res., 105, 
15959 

Belcher, J. W., & Davis, L. 1971, J. Geophys. Res., 76, 3534 

Belcher, J. W., Davis, L., & Smith, E. J. 1969, J. Geophys. Res., 74, 2302 

Breech, B., Matthaeus, W. H., Minnie, J., Bieber, J. W., Oughton, S., Smith, 

C. W., & Isenberg, P. A. 2008, J. Geophys. Res., 113, A08105 

Bruno, R., & Carbone, V. 2005, Living Rev. Sol. Phys., 2, 4 

Chandran, B. D. G. 2008, ApJ, 685, 646 

Chandran, B. D. G., Quataert, E., Howes, G. G, Hollweg, J. V., & Dorland, 

W. 2009, ApJ, 701,652 

Coleman, P. J. 1968, ApJ, 153, 371639, 1177 

Cranmer, S. R., Matthaeus, W. H., Breech, B. A., & Kasper, J. C. 2009, ApJ, 
702, 1604 

Cranmer, S. R. & van Ballegooijen, A. A. 2005, ApJS, 156, 265 (CvB05) 
Cranmer, S. R., van Ballegooijen, A. A., & Edgar, R. J. 2007, ApJS, 171, 520 
Del Zanna, L., Velli, M., & Londrillo, P. 2001, A&A, 367, 705 
Dmitruk, P., Milano, L. J., & Matthaeus, W. H. 2001, ApJ, 548, 482 
Dmitruk, P., Matthaeus, W. H.. Milano, L. J.. Oughton, S., Zank, G. P., & 
Mullan, D. J. 2002, ApJ, 575, 571 ("D02") 
Dmitruk, P., & Matthaeus, W. H. 2003, ApJ, 597, 1097 
Esser, R., & Sasselov, D. 1999, ApJ, 521, L145 

Feldman, W. C, Habbal, S. R., Hoogeveen, G., & Wang, Y.-M. 1997, J. 
Geophys. Res., 102, 26905 

Galeev, A. A., & Oraevskii, V. N. 1963, Sov. Phys.— Dokl., 7, 988 
Goldstein, M. L., Roberts, D. A., & Matthaeus, W. H. 1995, ARAA, 33, 283 
Grappin, R., & Velli, M. 1991, Ann. Geophys., 9, 416 



Grappin, R., Velli, M., & Mangeney, A. 1993, Phys. Rev. Lett., 70, 2190 

Heinemann, M.. & Olbert, S. 1980, J. Geophys. Res., 85, 1311 

Hollweg, J.V. (1983), Coronal heating by waves, in Solar Wind Five, edited 

by M. Neugebauer, p. 3, NASA Conf. Publ. 2280, Washington DC. 

Hollweg, J. V. 1986, J. Geophys. Res., 91, 4111 

Hollweg, J. V., & Isenberg, P. A. 2002, J. Geophys. Res., 107, 1 

Hollweg, J.V, & Isenberg, P. A. 2007, J. Geophys. Res., 112. A08102 

Hollweg, J. V, & W. Johnson 1988, J. Geophys. Res., 93, 9547 

Kadomtsev, B. B. & Pogutse, O. P. 1974, Sov. JETP, 38, 283 

Marsch, E., & Tu, C.-Y. 1996, Ann. Geophys. 14, 259 

Matthaeus, W. H., Zank, G. P., Oughton, S. , Mullan, D. J., & Dmitruk, P. 
1999, ApJ, 523, L93 

Roberts, D. A., Goldstein, M. L., Matthaeus, W. H., & Ghosh, S. 1992, J. 
Geophys. Res., 97, 17115 

Schekochihin, A. A., Cowley, S. C, Dorland, W., Hammett, G. W., Howes, 

G. G., Quataert, E., & Tatsuno, T. 2009, ApJS, 182, 310 

Spruit, H. 1981, "Magnetic flux tubes," in The Sun as a Star, ed. S. Jordan, p. 

385, NASA SP-450, Washington DC. 

Strauss, H. R. 1976, Phys. Fluids, 19, 134 

Tu, C.-Y. 1987, Solar Phys., 109, 149 

Tu, C.-Y. 1988, J. Geophys. Res., 93, 7 

Tu, C.-Y, & Marsch, E. 1995, Sp. Sci. Rev, 73, 1 

Tu, C.-Y, Pu, Z.-Y, & Wei, F.-S. 1984, J. Geophys. Res., 89, 9695 

Velli, M. 1993, Astron. Astrophys., 270, 304 

Velli, M., Grappin, R., & Mangeney, A. 1989, Phys. Rev. Lett., 63, 1807 
Velli, M., Grappin, R., & Mangeney, A. 1990, Comp. Phys. Comm., 59, 153 
Verdini, A., & Velli, M. 2007, ApJ, 662, 669 (VV07) 
Verdini, A., Velli, M., & Buchlin, E., 2009, ApJL, 700, 39 
Vifias, A. E, & Goldstein, M. L. 1991, J. Plasma Phys., 46, 129 
Zank, G. P. & Matthaeus, W. H. 1992, J. Plasm. Phys., 48, 85 



